Simulations of Magnetic Field Generation in Unmagnetized Plasmas via Beat Wave 

Current Drive 
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This work describes the scientific basis and associated simulation results for the magnetization of 
an unmagnetized plasma via beat wave current drive. Two-dimensional electromagnetic particle-in- 
cell simulations have been performed for a variety of angles between the injected waves to demon- 
strate beat wave generation in agreement with theoretical predictions of the beat-wave wave vector 
and saturation time, revealing new 2D effects. The simulations clearly demonstrate electron ac- 
celeration by the beat waves and resultant current drive and magnetic field generation. The basic 
process depends entirely on the angle between the parent waves and the ratio of the beat-wave phase 
velocity to the electron thermal velocity. The wave to magnetic energy conversion efficiency of the 
cases examined is as high as 0.2%. The technique could enable novel plasma experiments in which 
the use of magnetic coils is infeasible. 
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The nonlinear mixing of electromagnetic (EM) waves 
in plasmas has been applied to heating, diagnostics, par- 
ticle acceleration, ionospheric plasma modification, and 
current drive This paper describes the application 
of beat wave current drive [2|-|4j to magnetize an initially 
unmagnetized plasma. The greatest emphasis of prior 
beat wave studies was on accelerating electrons to rela- 
tivistically high energies and not on the most effective 
generation of current for producing significant magnetic 
field in an ambient plasma. Potential applications of beat 
wave magnetic field generation include controlled fusion 
concepts such as magneto- inertial fusion (MIF) [5H2J , and 
laboratory experiments on plasma phenomena of astro- 
physical interest, e.g., For lack of immediate appli- 
cations, however, the full complexity of significant cur- 
rent and field generation from beat wave current drive 
has not been addressed. While early analyses of beat 
wave current drive were essentially ID Q , more realistic 
modeling is now required with actual exploratory exper- 
iments underway [9J. For significant current drive, the 
beat wave phase velocity and electron thermal velocity 
must be comparable in order to accelerate a useful num- 
ber of electrons, and thus both collisionless (Vlasov) and 
collisional modeling are required. The modeling must 
also be at least 2D in order to handle essential exper- 
imental input parameters, such as the shapes, widths, 
and injection angle of the overlapping EM waves. It is 
also required to obtain the resultant spatial distribution 
of extended return currents and the magnetic field. The 
2D modeling also allows more realistic testing of predic- 
tions from ID analysis Q. Using the Large Scale Plasma 
(Lsp) code [lOj . we have performed such a study using 
a particular EM wavelength and plasma density range of 
relevance to MIF. 

Remote magnetization has been accomplished previ- 
ously by shining a laser onto foil targets with wire loops 
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frequency waves is negligible and placement of the beat 
wave interaction region within the plasma can be precise; 

(2) beat waves are produced in a controllable direction 
depending on the angle 9 between injected waves; and 

(3) current drive can be accomplished for thermal plas- 
mas via control of the wave phase velocity which also 
depends on 9. For MIF energy applications requiring re- 
duction of cross-field energy transport, the Hall parame- 
ter (ratio of the electron cyclotron to momentum trans- 
fer frequency) uj c /is m > 1. This ratio can be achieved 
by seeding modest magnetic fields (B ~ 1 T) in lower 
density regions of imploding plasmas that get amplified 
geometrically to ~ 100 T at maximum compression Q. 

Figure [T] shows the setup and representative simu- 
lation results using two CO2 lasers, with intensity of 
3 x 10 12 W/cm 2 and 9 = 90°, impinging on the center 
of a circular plasma with non-uniform density peaking 
at 3 x 10 16 cm~ 3 and initial temperature To = 50 eV. 
The orientation of the lasers and the subsequent beat 
wave are shown in Fig. 1(a), with the beat wave vec- 
tor being the difference between the two laser wave vec- 
tors, i.e., kb w = k\ — k,2- The upward- and rightward- 
propagating lasers have 10.4 and 10.8-/im wavelengths, 
respectively. The combined laser electric field pattern, 
shown in Fig. 1(b), drives a plasma beat wave propagat- 
ing to the upper left, as shown in Fig. 1(c). The beat 
wave modulates the plasma density by roughly ±10% 
with 7-/mi wavelength. The plasma beat wave acceler- 
ates a portion of the plasma electrons to as high as 500 eV 
in the kb w direction, driving a current and magnetic field 
structure shown in Fig. 1(d). 

We now discuss the theory of beat wave current drive 
as exhibited in Fig. 1. Two waves, with frequencies wi, 
u>2 3> the electron plasma frequency ui pe , are launched 
into a plasma to generate a beat wave with frequency 
13L 14]. Resonant interaction between the beat 



Beat wave magnetization, however, offers several 
potential advantages: (1) refraction of the injected high 



wave and the electron population is exploited to accel- 
erate electrons via Landau damping and to drive cur- 
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FIG. 1. Beat wave magnetic field generation using two lasers 
of 3 x 10 12 W/cm 2 intensity and 9 = 90° in a 50-eV plasma 
with 3 x 10 16 -cm -3 peak density, (a) Laser injection and beat 
wave geometry; (b) laser electric fields 20 ps after injection; 
(c) modulated plasma density (at 20 ps); and (d) B contours 
and electron mean velocity vectors (at 100 ps). 



rent and generate magnetic fields. Consider a beat 
wave generated from the ponderomotive force {E x B) 
of two intersecting injected waves with electric field E 
polarization into the page and wave vector as shown in 
Fig. 1(a). The angle a between kf, w and k\ is tan(a) = 
[&2 sin(0)]/[fci — fc 2 cos(f?)]. If we include the ponderomo- 
tive force in the cold electron fluid equation of motion, 
the magnitude of the wave vector is 



kbw — 



k\ — k 2 cos(0) 
cos(a) 



(1) 



Furthermore, the rate of increase in the beat wave am- 
plitude is [H[ 



At) 



kbw 

4w. 
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where ej = eEi/mcuji, and uii is the individual laser elec- 
tron oscillation frequency. For the case that the beat 
wave phase velocity v p h = u> pe /khw <C c, the wave will 
break for Akb w = 1 and saturate [14[. Thus, we can esti- 
mate the saturation time of the beat wave from Eqs. ([!]) 
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Although the beat wave growth rate scales linearly with 
laser intensity [ex (Ii^) 1 ^ 2 ], the wave saturated ampli- 
tude is a weak function of laser intensity. For cases pre- 
sented in this paper, r sat ~ 1 ps, providing a lower bound 
on the simulation duration. 

As shown in Fig. 1, we model the entire beat wave in- 
teraction and current drive using Lsp, a statc-of-thc-art, 
parallel, multi-dimensional particle-in-cell (PIC) code. 
The simulations resolve the smallest relevant EM wave- 
lengths (~ 5 /im) and highest frequencies (~ 8 THz). 
The > 0.1-ns duration and order mm scale length for 
the interaction demand 10 5 time steps and > 10 6 cells 
in 2D. Hundred of particles per cell are required to ad- 
equately resolve the electron energy distribution in the 
high energy tail and provide a sufficiently small noise 
level to discern the beat wave dynamics. For these simu- 
lations, Lsp solves the relativistic Maxwell-Lorentz equa- 
tions with inter-particle collisions. Particle scattering is 
treated with complete generality via a binary Coulombic 
interaction algorithm [IcJIlE}- We use an explicit, energy- 
conserving particle advance that sums particle currents 
such that charge is conserved [la ]. 

Our simulation domain is a 1-mm by 1-mm square 
enclosing a 0.95-mm diameter plasma cylinder (with 
vacuum between plasma edge and domain boundary) 
as shown in Fig. [TJc) with density range from 1-3 x 
10 16 cm -3 . From the boundaries, we inject 100-/im 
transverse-extent lasers into the plasma. Intensities ex- 
amined here are relevant for generating seed magnetic 
fields of order 5 kG. The resonant plasma density for the 
case of 10.4 and 10.8 /xm wavelengths is 1.5 x 10 16 cm" 3 , 
half the peak density. We found that this difference is 
not significant to the resulting beat wave production al- 
though plasma heating is more efficient at the resonance. 
The lasers reach 10 12 -10 13 W/cm 2 peak intensity in 10 ps 
with 9 = 20-180°. A range of beat wave characteristics 
can be obtained by varying 9. Much past theoretical 
work focused on co-linear (9 = 0) injection for its rele- 
vance to relativistic electron acceleration. In that case, 
Vph approaches c, which is suitable for a high energy ac- 
celerator, but does not easily couple to a thermal plasma 
electron distribution as required for current drive. There- 
fore, we focus most of this paper on larger 9. 

We have verified the simulations against basic theory. 
In all but the 9 = case, the beat wave amplitude sat- 
urates in several ps, so we present the saturated wave 
characteristics after 20 ps in several simulations. The 
beat wave characteristics including wavelength and direc- 
tion are nicely exhibited by the electron density contours 
in Fig. [2] [and Fig. 1(c) for 9 = 90°]. The simulations 
agree closely with the theory for both a and kbw- The 
smallest kbw is 1/3 that of the lasers for the 9 = 20° 



case, and the largest, roughly twice that of the lasers, is 
found for the 180° case. The beat wave saturated ampli- 
tude, as measured by the wave displacement A or elec- 
tron density modulation, is predicted to scale inversely 
with kb w . This is not obvious in the simulations largely 
due to the many other competing nonlinear interactions 
including temperature dispersion, collisionality, and cur- 
rent drive. The beat wave electron density modulations 
are all roughly ±10%. Peak beat wave electric field mag- 
nitudes Eb w range from 100-1000 kV/cm with only a 
weak dependence on intensity. 

We find that current drive is higher for larger angles 
(Fig. [2|. For To = 50 eV, line currents of 8 kA/cm pro- 
duce B ~ 2 kG. The 9 = 20° and 47° simulations show 
an order of magnitude smaller fields due to the fast v p h 
of the beat waves. In all cases, the electron current is 
driven in the direction of kb w . The mean electron ve- 
locity vectors in Fig. [2] reveal a closed current path sur- 
rounding magnetic islands for the 9 = 90-180° simula- 
tions. The lower-0 simulations have a less obvious re- 
turn current path due to the finite volume of the simula- 
tion plasma. Higher energy electrons carrying the current 
quickly reach the plasma-vacuum interface. 

Although the beat wave saturation amplitude is only 
a function of 9, the beat wave current drive is optimized 
by controlling F = v p h/v te (where v te is the electron 
thermal velocity), which varies in time by plasma heat- 
ing observed at higher laser intensity. In Fig. [3ja), we 
show peak B(t) for several 9 = 90° simulations with 
3x 10 12 W/cm 2 laser intensity and F varying from 1.3-11 
(To = 200-25 eV). The simulations show a rapid rise in 
B within 10 ps of peak laser intensity at the interaction 
region. The field saturation time is consistent with the 
theoretical value of T sa t ■ An optimal value for F is in the 
range 1.9-2.7 for which S ~ 1 kG. 

The beat wave current drive efficiency is [3] 

[E ph /\Q keV) 1 - 5 
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2(L int /10 m) ne/10 1 



-(A/W), (4) 



where Li nt is the laser interaction length (roughly 5 cm 
in our case). A conservative estimate of the action trans- 
fers i?i and i? 2 of 0.1 is taken from Fig. 1 of Ref. 3, 
and E p h is the electron energy associated with v p h- Us- 
ing Eq. ((U) to estimate peak B, we compare the theory 
with simulation results as a function of laser intensity for 
9 > 90°. As shown in Fig. [3jb), the ID theory and 2D 
simulations follow a similar linear dependence of B on 
intensity, although there is some variation due to vary- 
ing F. We also calculate the total conversion efficiency 
of laser to magnetic field energy which we define as the 
ratio of magnetic field energy at 100 ps to the total in- 
jected laser energy Ei n j through 100 ps. Keeping fixed 
9 = 90° and 3 x 10 12 W/cm 2 intensity (E inj = 0.015 J), 
we find a maximum in the energy efficiency of 0.03% near 
F = 2.2 with r] > 0.001% for F = 1.9-2.7, as shown in 
the inset of Fig. IS^a). The efficiency falls off steeply for 
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FIG. 2. Beat wave modulated electron density contours (left 
column) at t = 20 ps for (a) 20°, (c) 47°, (e) 135°, and 
(g) 180°; magnetic field contours and electron mean veloc- 
ity vectors (right column) for (b) 20° at 50 ps, (d) 47° at 50 
ps, (f) 135° at 100 ps, and (h) 180° at 100 ps. This is for 
laser intensities of 3 x 10 12 W/cm 2 and 50-eV plasma. 



F > 3. The highest efficiency of 0.2% seen in all our sim- 
ulations thus far correspond to 9 = 180° and 10 13 W/cm 2 
(Emj — 0.045 J). The efficiency can be optimized via F 
and maximizing the overlap of injected beams. 

Examining more closely the mechanism of current 
drive in our nominal 3 x 10 12 W/cm 2 intensity and 
6 = 90° simulation where v p h — 0.027c, we find that 
beat waves with peak Eb w > 400 kV/cm can trap elec- 
trons moving in the same direction with sufficient initial 
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FIG. 3. Peak magnetic field driven by beat waves versus 
(a) time at 3 x 10 12 W/cm 2 and 6 = 90° for different values 
of F = v p h/vte, and (b) laser intensity for 9 > 90° including 
theoretical values calculated from Eq. Q and a linear fit to 
the simulation data. The inset in (a) shows the efficiency of 
magnetic energy generation versus F at 100 ps. 




08.0 0.03 0-04 0.05 0.06 0.07 0.03 0.04 0.05 0.06 0.07 





0.03 0.04 0.05 0.06 0.07 0.03 0.04 0.05 0.06 0.07 

Z (cm) Z (cm) 

(c) (d) 



FIG. 4. Electric field magnitude at t — 50 ps for 3 x 
10 12 W/cm 2 laser intensity, 6 — 90°, and varying To: 
(a) 10 eV (F = 11), (b) 25 eV {F = 3.8), (c) 50 eV (F = 2.7), 
and (d) 75 eV (F = 2.2). 



speed, such that v > v ph - [AE bw / (mc 2 k bw )] 1 / 2 s» 0.008c 
(16 eV), and accelerate them to higher energy. This ve- 
locity is determined from the total depth of the poten- 
tial well of the wave in the frame of v p h- We injected 
test electrons in the center of the laser interaction region 
(X = Z = 0.05 cm) with (0.0014-0.02)c velocity in the 
direction of the beat wave and tracked their velocity ( V) 
and position (X') in that direction. Electrons with ini- 
tial energy < 18 eV exhibited oscillatory behavior with 
peak energy increasing to 50 eV, but little net motion. 
Most higher energy electrons are trapped, riding up and 
down on the beat wave fronts but moving with the beat 
wave. The electrons accelerate to V « 0.045c (500 eV), 
then down to their initial values for several cycles until 
they are scattered or begin to leave the region of laser 
interaction at X' > 0.01 cm. At this time, the elec- 
trons can retain some fraction of their peak accelerated 
energy. These strongly-directed electrons drive current 
until they are scattered or reach the edge of the plasma. 
Some of these electrons end up carrying return current as 
they cycle back. The detailed relationship between beat 
wave fields, electron acceleration, and magnetic field ex- 
tent will be explored further in a forthcoming paper. 

The local E\, w , in regions of highest electron acceler- 
ation and at optimal F, decays due to nonlinear Lan- 
dau damping. The Landau damping rate is significant as 
kbw^D 1- For the F = 2.2 simulation, kbw^D = 0.3. 
The damping effect on beat wave fields is shown in Fig. [4] 
for varying To. For this particular beat wave excitation 
at a given phase velocity, at low temperature [Fig. 4(a)] 
when Landau damping is not strong, the structure re- 
flects quite well the beat wave excitation with negligi- 
ble loss. With higher temperature [Fig. 4(d)], significant 



damping is now evident in the decay of the beat wave 
excitation magnitude along the beat wave wavevector di- 
rection (at > 135° to the Z axis). Indeed at 200 eV, 
the beat wave is strongly damped everywhere relative to 
the 10-eV simulation. It is unclear whether the presence 
of the magnetic fields themselves, producing cyclotron 
radii « 0.02 cm, also contributes to the decay of the beat 
wave. The regions of stronger beat wave damping also 
have the highest magnetic field. The 200-eV simulation 
had 1/5 the peak magnetic field but faster damping of 
the wave than the 75-eV case. The Landau damping of 
the beat waves permits the accelerated electrons to es- 
cape the wave structure and drive current. Thus, the 
strong wave and current regions remain distinct. 

In this paper, we have described the scientific basis for 
magnetization of an unmagnetized plasma via beat wave 
current drive. Specifically, we performed 2D PIC simu- 
lations to explore the scaling of beat wave production, 
current drive, and B generation in a non-uniform den- 
sity plasma using injected waves near 10-/j,m wavelength, 
corresponding to CO2 lasers. Future work should include 
3D simulations and exploration of lower intensity lasers 
(10 8 W/cm 2 ) for modeling of near term experiments 0[. 
The 2D simulations to 100 ps presented here made use 
of 64 processors run for roughly 72 hours. Simulations in 
3D will be 100's of times larger in cell number and require 
1000's of processors on a massively parallel computer. 

We acknowledge excellent code support from R. E. 
Clark. This work was supported by the Office of Fusion 
Energy Sciences of the U.S. Department of Energy. 
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